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The idea of writing a table of probabilistic data for a quantum or classical system, 
and of decomposing this table in a compact way, leads to a shortcut for Hardy's 
formalism, and gives new perspectives on foundational issues. 
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Dcnkbarcs, Mcnschcn-Sichtbarcs, Menschcn- 
Ftihlbares! Eure eigncn Sinnc sollt ihr zu Endc 
dcnkcn! 

Nietzsche 



1 Introduction 

Suppose that we have written, in a sort of table, the statistical data collected 
from a group of experiments — the nature of which can be classical, quantum, 
or something else. Suppose that we also want to store this table's data in a 
compact way. How could we proceed? 

In this paper it is shown that, given the situation described above, when 
we try to store or organise the table's data in a more compact way we find that 
real vectors can be associated to preparations and results, in a way which, for 
quantum mechanical phenomena, is essentially the same as Hardy's represen- 
tation of 'states' and 'measurement outcomes' This curious fact may offer 
new points of view for looking at some of the current 'foundational' issues in 
quantum mechanics. 

The ideas here presented are a summary of those developed in Ref. |2|, to 
which the Reader is referred for further details. The emphasis in this paper 
is on the main idea of a 'table decomposition', and on the implications of the 
latter for various topics discussed in this Conference. a 



a Moreover, many analogies, discussed in Ref.0, are to b e found between t his work and those 
of Mielnik QUI 00 Foulis and Randall et al. B 0, E3, E3 , Barnum 0, and others. 
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2 Probability tables 



Imagine that we are in a laboratory, performing experiments of various kinds 
to study some interesting phenomena; the purpose of the experiments is to 
statistically study the correlations among different kinds of these phenomena. 
In general, we try to reproduce a given phenomenon — either by controllably 
preparing it at will, or simply by waiting for its occurrence — , to observe which 
concomitant phenomena, or results, occur. 

Some experiments present common features: for example, part of the 
preparation can be the same for some of them. We separate ideally each 
experiment into a preparation and an intervention; the latter also delimits the 
kind of results which can be obtained, which implies that if we are told a result, 
we know which intervention was made. We then consider a set of preparations 
and a set of interventions, with the clause that sensible experiments may be 
made by combining each of the preparations with each of the interventions 
(preparations or interventions which do not satisfy this condition are set aside 
for the moment ). b 

Thus, suppose that we have M different preparations {Si, . . . , Sm}, and 
a given number of possible interventions {Mi, M 2 , ■ ■ ■ , Mk, ■ ■ ■}, each with a 
different number of results {R\, R 2 , ■ ■ ■ , Ri M } (mutually exclusive and 
exhaustive ), where the number lM k depends on the particular intervention 
Mfe. The total number of results, counted from all interventions, is L. 

Through repetitions of the experiments, or through theoretical assump- 
tions, or just by analogy with other experiments which we have already seen 
and which we judge similar to those that are now under study, we can write 
down a table p with the probabilities that we assign to every result, for every 
intervention and preparation. The table may look like the following: 



Note that a preparation does not need to temporally precede an intervention; indeed, these 
two notions arc meant to have here only a logical, not temporal, meaning. For example, 
in quantum-mechanical experiments with post-selection, the preparation is effectively com- 
pleted after the intervention is made! 

c This can always be achieved by grouping in suitable ways the results, and adding if necessary 
the result "other" . 
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P64 
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Rl 
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PL2 


PL3 


Phi 


PLM 



The table, which may be called a 'probability table', has a column for each 
preparation and a group of rows for each intervention, and these rows are the 
possible results of the intervention. The table entry p^ is the probability of 
obtaining the result Ri for the intervention and the preparation Sj. For 
example, the entry (4, 3) is the probability p4 3 the we assign to obtaining the 
result /?4, among the possible results {R 3 , R4, R5}, for the intervention M2 
and the preparation S3 . Preparations and results can be listed and rearranged 
in any desired way in the table. Such a table would very likely have a large 
number of rows and columns, i.e., the numbers L and M are likely to be very 
large. 

Now, suppose that we seek a more compact way to write down and store 
the probability data collected in the table p. We note that the table is really 
just an L x M rectangular matrix, and as such it has a rank K, viz., the 
minimum number of linearly independent rows or columns: 



K = f rankp ^ min{i, M}. 



(1) 



It follows from linear algebra that p can be written as the product of an L x K 
matrix t and a K x M matrix u: d 



p = tu, 



(2) 



or 



Pll ■ 


■ Pli ■ 


■ PIM \ 




('I 


Pil ■ 


■ Pij ■ 


• PiM 






PLl ■ 


■ PLj ■ 


■ PLM / 




u 



(3) 



In the last equation, the matrix t has been written as a block of row vectors 
rj, and the matrix u as a block of column vectors Sj. In this decomposition, 



d This is equivalent to the fact that a linear map p: 



WL L of rank K d = dimp(lR M ) 



can be obtained as the composition p ■■ 
injective map t: p( 



t o u of a surjective map u: R M — ► p(R M ) and an 
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the element p^ of p is then given by the matrix product of the row vector rj 
with the column vector s j : 

Pij =rjsj =ri-Sj, (4) 



where, in the last expression, and Sj are considered as vectors in R , so 
that the matrix product is equivalent to the scalar product. It will be shown 
in a moment that the decomposition is always effective in reducing the number 
of data of the table. 

The result is that we can associate vectors {sj} and {r,} in H K , for some 
K, to the preparations and the intervention results for the table, and the 
relative probabilities are given by their scalar product: 



These vectors can be called preparation vectors and (intervention- Jresult vec- 
tors, or, in general, (table) vectors. 

It should be remarked immediately that we are not postulating any kind 
of physical property in Eq. (JSJ); even less have we found one. We have only 
decided to represent and store a collection of numbers (which do have physical 
significance) in an alternative way. Note, in particular, that the numerical 
values of the vectors {sj} and {r^}, as well as their dimension K, depend on 
the whole collection of probabilities {pij}'- if one of these is changed, then K 
and all the vectors will in general change. 

The meaning of the representation (J5J shall be discussed in a moment, but 
let us study the decomposition in more detail first, in order to have a clearer 
idea of how the table vectors originate. 

The matrices t and u are not uniquely determined from the decomposi- 
tion J2J, so that there is some freedom in choosing their form. The fact that 
rankp = K, implies that there exists a square K x K submatrix a, obtained 
from p by suppressing (L — K) rows and (M — K) columns, such that det a ^ 0. 
It is always possible to rearrange the rows and the columns of the table p so 
that this non-singular submatrix is the one formed by the first K rows and 
K columns. After this rearrangement, p can be written in the following block 
form: 



where b, c, and d are of order Kx(M-K), (L-K) x K, and (L-K)x(M-K) 
respectively. 

By writing also the matrices t and u in block form 



Pij 



(5) 




with det o^O, 



(6) 




(7) 
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where the orders of v, w, x, and y arc K x K, (L — K) x K, K x K, and 
K x (M — K) respectively, we can rewrite the decomposition equation (J2j) as 



(: 5) -(:)(■»> ■ 



(8) 



or 



a = vx, b = vy, c = wx, d = wy. 



(9) 



which has the solution 



detx^=0, y = xa 1 6, v = ax , w = cx 



l 



(10) 



which in terms of the result and preparation vectors is 



(si ... s K ) = x, (s K+1 . . . s M ) = xa 1 6, 



(11a) 





(lib) 



The square matrix x = (si . . . sk) is undetermined except for the condition of 
being non-singular; this corresponds to the freedom of choosing K basis vectors 
in R x as those associated to the first K preparations f Si, ... , Sk, which can 
then be called basis preparation-vectors. Note that for a given probability 
table p the set of basis preparation-vectors is in general non-unique, because 
p possesses in general many submatrices like a with rank K . 

We can now check whether the above decomposition effectively reduces the 
number of data to be stored. The table p has L x M entries. The matrices x, 
y, v, w (or equivalently the collection of vectors {sj} and {t^}) have in total 
K x (L + M) entries; however, we are free to choose a 'canonical' form for the 
non-singular K x K matrix x (e.g., the identity matrix) once and for all, for 
all probability tables (this corresponds to a standard choice of basis vectors 
in R K ). This implies that we are indeed left with only K x [L + M) — K 2 
numbers. It is then simple to verify that, from K ^ min{L, M}, we have 



Thus the decomposition of the probability table into two sets of vectors is 
indeed a more compact way to write down and store the probability data. 

e Note that d = car b, i.e., the submatrix d of p is completely determined by the other 
submatrices a, 6, and c, because rankp = K. 

f There is the alternative option of choosing the vectors associated to the first K results; this 
corresponds to solving Eq. J2J in terms of v. 



K x (L + M) - K 2 ^ Lx M. 



(12) 
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2.1 A numerical and graphical example 

Imagine that, with the various apparatus in our laboratory, we can make seven 
different preparations {Si, . . . , S7} and perform three different interventions 
{Mi, M 2 , M 3 }, each having two results. The probabilities that we assign to 
the various results are given in the following table: 





Si S2 S3 S4 S5 S6 S7 




1 i n 1 3 I 3 

1 2 U 2 4 2 4 

I l I I I I 




T ? 1 1 3 ? r— 

2 1 2 u 4 2 2 

1 - 1 i - i 

2 2 4 2 2 


M > s 


1111111 





The matrix corresponding to this probability table has rank X = 3, and 
the first 3x3 submatrix a is indeed already non-singular, without the need 
of rearranging the rows or columns of the original table. Proceeding as in the 
decomposition equations JHJ and (|10|l with the following choice for the matrix 
x: 

'II 1 

10-lJ, (13) 
v 1 

we finally obtain the following preparation and result vectors: 



S\ 


= (1,1,0), 


S-2 


= (1,0,1), 


S3 


= (1,-1,0), 


S 4 


= (1,0,-1), 


S5 


= fl 1 1) 


s 6 


= (1,0,0), 


st 


= (1.|,Q). 



n = (3,3.0), 

r-2 = (i-|,0), 
r 3 = (!,0,i), 
r 4 = (1,0,-1), 
r 6 = (1,0,0), 
r 6 = = (0,0,0), 



(14) 



represented as points of R 3 in Fig. 

We notice that all preparation vectors lie on the same plane (x = 1); this 
is a general property of any probability table, which derives from the fact 
that the results of an intervention are mutually exclusive and exhaustive; this 
also implies that an intervention's result-vectors always sum up to the same 
vector |2|: in this case, r\ + r 2 = r 3 + r 4 = r 5 + r 6 = (1,0, 0). 
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3 Relation to quantum mechanics 



Probability tables like the one illustrated above can be made, in particular, for 
phenomena concerning classical or quantum systems. Indeed, the resulting as- 
sociation of vectors to preparations and results is analogous to that introduced 
by Hardy 0, through a different line of reasoning, for 'states' and 'measure- 
ment outcomes' of classical and quantum systems. 

For example, imagine how a table for a two-level quantum system would 
appear; consider for concreteness the polarisation of a single photon. We can 
prepare the photon with different polarisation directions and with different 
degrees of polarisation; hence the table has in the limit a continuum of columns, 
one for each preparation. Analogously, we can perform some interventions on 
the photon by placing various polarisation filters on its way, controlling if it is 
absorbed or not; in the limit, the table has also a continuum of rows, one for 
each intervention result. Some entries of the table would look like the following 
(the meaning of the symbols should be self-evident): 
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30° ■ ■ ■ ~>45° 






^45° 045° 
K abs 


... \ ... 1 

... | ... 






p60° 

A/f "out 
'^60° p60 ° 

"abs 


\ ... 0.933 
§ ... 0.067 







Such a table has, notwithstanding the limiting infinite size, rank K = 4, 
and so we can associate to every preparation and to every result a 4-dimensional 
vector. The preparation vectors, however, lie on a 3-dimensional (afhne) hyper- 
plane, as in the numerical example previously discussed. It can be shown 
indeed 0, 0] that the resulting set of preparation vectors is equivalent to the 
standard Bloch-sphere for two- level quantum systems. 

This was just an example, but all quantum-mechanical concepts like den- 
sity matrix, positive-operator-valued measure, and completely positive map 
can be expressed in an equivalent table- vector formalism 00 (e.g., the rela- 
tion with the trace rule is quickly shown in the appendix). 

4 Discussion 

The following discussion concentrates on the possible relations among the ideas 
hitherto presented and ideas presented by other authors in this Conference. 
The Reader is referred to Ref. |2| for a more general and detailed discussion. 

4-1 On the probabilities in the table and 'quantum logics' 

In Sect. [21 we quickly introduced the preparations, interventions, and results 
{Sj}, {Mk}, {R-i}, and the probabilities {pij} which make up the table. Let 
us define them more clearly. 

Symbols like Si and Mk represent propositions which together describe an 
actual, well-defined procedure to set up an experiment, e.g. 5; ='The diode 
laser is put on the table is such and such position, with a beam splitter in front 
of it in such and such position. . . ' etc., 8 and Mk = 'Such and such vertical filter 
is placed in such and such place, and the detector is placed behind it. . . ' etc. 
The separation of the experiment's description into the two propositions is not 
unique, and indeed more kinds of separations can be considered. 11 

g Anothcr example: 'Wait until such and such happens, then. . . ' etc. 

h In Ref. 0, this fact is used as the starting point to define the concept of transformation. 
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A symbol like Rj represents a proposition which describes the results of 
an experiment, e.g. Rj = 'The detector does not click'; it is then clear that it 
depends on the particular experiment being performed. 

The probability p^ is consequently defined as 



where Q is a proposition representing the rest of the experimental details and 
our prior knowledge. 1 We are thus considering probability theory as extended 
logic an approach which will prove to be, in the following, 

powerful, flexible, and intuitive at once. 

Usually, many repetitions of the same experiment are made and the relative 
frequencies of the different results of the intervention are observed. In this case, 
given a judgement of total exchangeability of the experiment 's repetitions, the 
probability is practically equal to the observed frequency of the result, thanks 
to de Finetti's representation theorem But the probability can also be 

assigned on grounds of similarity with other experiments, or just by theoretical 
assumptions. 

Note that the table is silent with regard to the probabilities for the prepa- 
rations or the interventions, P(Sj\ Q) and P(M k \ Q)- In a given 'situation' Q, 
if we decide which preparation and intervention to perform, say S5 and M7, 
then these probabilities are P(5;| Q) — 6^5 and P(Mk\ Q) — 6fe,7 of course 
(which amounts to saying that "we know what we're doing"). If it is someone 
else who is deciding the particular preparation and intervention, then the prob- 
abilities must be assigned in some other way, e.g. by asking "which preparation 
are you making?" , or by using some other knowledge, and they will in general 
differ from and 1. The difference between these probabilities and those of 
Eq. (|15|l is somehow analogous to the difference between initial conditions and 
equations of motion in classical mechanics: the theory concerns only the latter, 
while the former has to be specified on a case-by-case basis. 

The probabilities of the results of a given intervention and a given prepa- 
ration form a probability distribution, because the results were arranged so as 
to be mutually exclusive and exhaustive. This implies that, for two results R k 
and R' k of a given intervention M k and for a given preparation Sj, we also have 
the trivial identity 



P(R k V R' k \M k A Sj A Q) = P(R k \ M k A Sj A Q) + P(R' k \ M k A Sj A Q), 



Pi j = P(Ri \ M ki A Sj A Q), 



(15) 



{r k +r' k )-Sj. 



(16) 



'A 'subjective Bayesian' can see Q as representing the 'agent'. 
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But with probability theory as logic we can also evaluate, for a given prepa- 
ration Sj, the disjoint probability for the results R and R" of two different 
interventions M' and M", just using the product and sum rules: 

P[R' V R"\ (M' V M") A Sj A (?] = P[{R' A M') V (R" A M")\ Sj A (?], 

= P(R'\ M' A Sj A (?) x P(M'\ Sj A (?) + 
P(R"\ M" A Sj A (?) x P(M"| S 3 A (?), 

= [r'P(M'| Sj A (?) + r"P(M"| Sj A (?)] • «,-, 

(17) 

where it is assumed that P(M' V M"\ Sj -A (?) = 1, i.e., we are sure that one or 
the other intervention was performed. 

The content of the formula above is intuitive: the occurrence of the result 
R implies that the intervention M was performed, and so analogously for 
the result R" and the intervention M" . Then the probability of getting the 
one or the other result depends in turn on the probability that the one or the 
other intervention was made, hence these probabilities appear in the last line 
of the above equation. However, as already said, the probabilities of these 
interventions are not contained in the table, but must be given on a case-by- 
case basis. 

As a result, the two disjoint probabilities in Eqs. (|16fl and l|17|) 'behave' dif- 
ferently, and the reason for this is intuitively clear. However, we can partially 
trace in this fact the source of much research and discussion on partially or- 
dered lattices and quantum logics for the set of intervention results |l£|,|lj), 20]. 
Roughly speaking, the point is that, for a classical system, there is the theoret- 
ical possibility of joining all possible interventions (measurements) in a single 
"total intervention" : the table associated to a classical system can then be 
considered as having only one intervention; thus one needs never consider the 
case of Eq. (I17JI . However, such "total intervention" is excluded in quantum 
mechanics, and one is thus forced to consider the case of Eq. (|17ll . 

From the point of view of probability theory as logic, instead, there is 
no need for non-Boolean structures thanks to the possibility of changing and 
adapting, by means of Bayes' theorem, the context (also called 'prior knowl- 
edge' [l(j or 'data' [ijj) of a probability, i.e., the proposition to the right of 
the conditional symbol 



Jin contrast, 'Kolmogorovian' probability, with its scarce flex ibil i ty w ith respect to contexts, 
deals with these features with difficulty. Loubenets' efforts |2 ll |22| are directed to amelio- 
rating this situation. 
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4-2 Unknown preparations and their 'tomography' 

On the other hand, we may have the following scenario: we have repeated 
instances of a given preparation, but we do not know which. By performing 
interventions on these instances and observing the results, we can estimate 
which preparation is being made. Introducing the proposition D representing 
the results thus obtained, we have 

P(S-\DAO)- P(D\S^Q)P(S ] \Q) , . 

This is a standard "inverse-inference" result of Bayesian analysis [3 Ch. 4]. 
The probability P(D\ Sj A Q) can be written in terms of scalar products of 
result and preparation vectors, but the probability distribution {P(Sj\Q)} 
depends on the prior knowledge that one has in each specific case. If we now 
want to predict which result will occur in a new intervention M^, we have the 
following probability: 

P(Ri\ M kt ADAQ) = Y, P(Ri\ M kl A Sj A Q) x P(Sj | D A Q), 

3 (19) 
= r i -Y d S j P{S j \DhQ), 

3 

i.e., we can effectively associate the vector 

«ncw = sj P(Sj \ D A Q) (20) 

3 

to the unknown preparation. 

A further kind of scenario is this: we have a brand-new kind of preparation, 
i.e., a new phenomenon, still untested. It has no place in our probability table; 
yet we think that we could reserve a new column to it without making substantial 
changes to our table's decomposition (by which we mean that the table's rank 
K would not change)^ Given this, in order to associate a vector to this new 
preparation we proceed as in the preceding scenario, performing interventions 
and observing results. The result is that Eqs. I|18|) and (|20|l also apply in this 
case. 

Note that there is no conflict between our talking about "unknow n prepa - 
rations" and Fuchs' criticism of the term 'unknown quantum state' |23u24l . l25| . 



k This condition seems to be somehow related to Fuchs' "accepting quantum mechanics" |23 
Sec. 6]. 
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A preparation, as we have seen, is a well-defined procedure (that can be shown 
or described to others, etc.) to set up a given physical situation; on the other 
hand, Fuchs' meaning of 'quantum state' is 'density matrix' [23L Eif. which 
corresponds (more or less) to the preparation vector instead. 1 Thus, consider 
the two sentences "It is unknown to me which kind of laser and of beam splitter 
are used in this experiment" and "I do not know with which probability the 
detector behind the vertical filter will click" : the latter sentence is nonsensical 
from a Bayesian point of view, because there are no 'unknown' degrees of be- 
lief; but the former sentence is unquestionably meaningful. Thus, even if the 
preparation is unknown to us, we can always associate a preparation vector to 
it. m 

4-3 On the content of the table vectors and 'quantum states of knowledge' 

Some remarks have already been made in Sect. after the derivation of the 
formula 

P(Ri\ Mki A Sj A Q) = tv sj, (ULis) 

with regard to the fact that the vectors {sj} and {rj} do not have any physical 
meaning separately: a given preparation vector Sj tells us nothing if we do not 
have an intervention-result vector r^; even less if we know nothing about the 
set of intervention results. It follows that the preparation vectors also lack 
any probabilistic meaning per se: they are not probabilities nor collections of 
probabilities — they are just mathematical objects which yield probabilities 
when combined in a given way with objects of similar kind. This, in particular, 
is true for density matrices, as we have seen that they are just a particular case 
of preparation vectors. It is slightly incorrect, as well, to say that probabilities 
are 'encoded' or 'contained' in the (quantum) preparation vectors: rather, they 
are parts of an encoding. 

From these considerations, the quantum state (the density matrix) appears 
to be part of a state of belief, and not the 'whole' state of belief jUH^]. Perhaps 
the point is that Caves, Fuchs, and Schack's notion of a 'quantum' state of belief 
implicitly assumes the existence and the particular structure of the whole set 

'Sadly, the Janus-faced term 'state' is sometimes meant as 'density matrix' and sometimes 
as a sort of 'physical state of affairs', though the two notions are, of course, quite distinct. 
m It would be interesting, indeed, to see more clearly the connexion between the "preparation 
tomography" illustrated above, and the 'q uantum Bayes rule ' and 'quantum de Finetti the- 
1 .ran ' of Caves, Fuchs, and Schack I2I |24 123. l27t I2I l2St l3fl ; it seems reasonable to expect 
that these 'quantum' analogues of Bayesian formulae should have a counterpart for tables 
concerning general systems, not only quantum-mechanical ones; the failure of the 'quantum 
de Finetti theorem' for quantum mechanics on real Hilbert spaces is quite interesting from 
this point of view. 
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of quantum positive-operator- valued measures (i.e., the interventions). 11 This 
is an important difference from a 'usual' degree of belief P(A\ X) which does 
not need to be combined with other mathematical objects to reveal its content. 



4-4 Possible applications of the table formalism 

It has already been remarked that the kind of vector representation arising from 
the table decomposition is essentially the same as Hardy's Q. The derivation 
presented here can be seen as a sort of shortcut for his derivation, but it implies 
something more. Hardy supposes that it is possible to represent a preparation 
by means of a i^-dimensional vector with K ^ L because most physical the- 
ories have some structure which relates different measured quantities; but the 
reasoning behind the decomposition of Sect.|21shows that this possibility exists 
even without a theory that describes the data (indeed, the question arises: has 
this possibility any physical meaning at all?) 

In any case, the idea of a 'probability table' and its decomposition has 
probably very little usefulness in experimental applications, but provides a 
very simple approach to study the mathematical and geometrical structures 
of classical and quantum theories, and offers a different way to look at their 
"foundational" and "interpretative" issues. 

This approach is even more general than other standard ones based, e.g., on 
C*-algebras, or even Jordan-Banach-algebras 31], or on convex state-spaces. 
Thus, with the idea of a 'probability table' we can very easily implement 'toy 
theories' or models like those of Spekkens and Kirkpatrick [33| (cf. also 
the issue raised by Terno [3), which can then be compared to classical or 
quantum mechanics using a unique, common formalism. 



"However, Caves, Fuchs, and Schack have the right to be not so pedantic about this point, 
because they assume at the start the absolute validity of quantum mechanics and of its 
mathematical structure (an assumption not made in the present work). 

°The first two formalisms are included as particular limiting cases of tables with an uncount- 
ably infinite number of columns and rows. In contrast to the convex-state-space approach, 
the probability-table idea does not require that the set of preparations be necessarily convex 
(however, the convexity usually appears in a natural way: see e.g. Hardy's discussion in 
Ref. Sect. 6.5). Moreover, in contrast to all three mentioned frameworks, it does not as- 
sume that the set of preparations (states) be necessarily the whole set of normalised positive 
linear functionals of the set of results (POVM elements), or vice versa for the convex-state- 
space framework. For example, (the convex hull of) the set of preparations for the table of 
Sect. I2T1 is only part of the set of normalised positive linear functionals of (the convex hull 
of) the set of results (the latter would be a square circumscribed on the given one). 
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A The trace rule 

It is shown that the 'scalar product formula', Eq. (JSJ), includes also the 'trace 
rule' of quantum mechanics (see also Ref. Q, Sect. 5). 

A preparation is usually associated in quantum mechanics to a density 
matrix pj, and a intervention result to a positive-operator-valued-measure ele- 
ment tli] both are Hermitian operators in a Hilbert space of dimension N. The 
probability of obtaining the result Tli for a given intervention on preparation 
Pj is given by the trace formula 

Pij = txtlipj. (21) 

The Hermitian operators form a linear space of real dimension K = N 2 \ one 
can choose K linearly independent Hermitian operators {Bk} as a basis for 
this linear space. These can also be chosen (basically by Gram-Schmidt or- 
thonormalisation) to satisfy 

tr B k Bi=b H . (22) 
Both pj and Ili can be written as a linear combination of the basis operators: 

K K 

P] = Si l B u = n k B k , (23) 

1=1 k=l 

where the coefficients s/ and ri k are real. Using Eqs. (|23|l and (|22|l the trace 
formula becomes 

K K 
Pij = tvflipj = ^2 r i ks 3 l tTB kBl = ^2 n l Sj l = Ti ■ Sj, (24) 
k,l=l 1 = 1 

where = (n 1 , ■ ■ ■ , ri K ) and Sj = (sj 1 , . . . , Sj K ) are vectors in R K , q.e.d. 
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